High pressure structures of disilane and their superconducting properties 
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A systematic ah initio search for low enthalpy phases of disilane (Si2H6) at high pressures was 
performed based on the minima hopping method. We found a novel metallic phase of disilane with 
Cmcm symmetry, which is enthalpically more favorable than the recently proposed structures of 
disilane up to 280 GPa, but revealing compositional instability below 190 GPa. The Cmcm phase 
has a moderate electron-phonon coupling yielding a superconducting transition temperature Tc of 
around 20 K at 100 GPa, decreasing to 13 K at 220 GPa. These values are an order of magnitude 
smaller than previously predicted Tc for disilane, and cast strong doubts on the possibility of high- Tc 
superconductivity in these systems as well as in other hydrogen-rich compounds under moderate 
pressure. 



Superconductivity in elemental hydrogen was pre- 
dicted by Ashcroft [l[ already in 1968. More recently, 
and with the use of novel theoretical techniques [2, 'H, 
the calculated Tc was estimated to be as high as 240 K at 
pressures of around 450 GPa However, the synthesis 
of metallic hydrogen has been found to be experimentally 
challenging, and even at extremely high pressures (below 
342 GPa) metallization has not yet been observed [HI, [6| . 
This is in agreement with theoretical calculations, that 
predict the metallic transition above 400 GPa — a pres- 
sure beyond the reach of current experimental capabili- 
ties. 

To circumvent this problem, it was recently suggested 
that metallization pressures could be achieved in hy- 
drogen rich materials where the hydrogen is chemically 
"pre-compressed" [7] . Several investigations of such com- 
pounds have appeared in the literature, primarily focus- 
ing on group IV-hydrides. Calculations on phases of 
highly compressed silane germane [14] and stan- 

nane ISl^llSl have shown the possibility of metallic phases 



with high Tc at moderate pressures. From the experimen- 
tal point of view, silane SiH4 has been reported to crys- 



tallize and attain metallicity above 50-60 GPa [17|, [18 1 
with a superconducting behavior. However, more recent 
studies ascribe the observed metallicity to the formation 
of metal hydrides [HI, and metallization of silane was 
found not to occur at least below 130 GPa (ioj . 

Another hydrogen rich compound of the same fam- 
ily is disilane Si2H6. This compound has attracted at- 
tention as a hydrogen rich material due to its exper- 
imental availability. Moreover, in a recent theoretical 
study, Jin et al. [21] performed random searches in or- 
der to find stable structures of disilane. They reported 
three different structures covering a pressure range from 
50 to 400 GPa. Crystallization of disilane into a metal- 
lic phase with a P-1 lattice was predicted to occur at 
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FIG. 1. (Color online) The crystal structures of (a) the Cmcm 
phase at 200 GPa, and (b) the P-1 phase at 300 GPa. The 
eigendisplacements which lead from the Cmcm structure to 
the Cmc2i structure are visualized by arrows in panel (c). 



135 GPa. The Tc of this phase was predicted to be 64.6 K 
at 175 GPa, and 80.1 K at 200 GPa. Beyond 175 GPa, the 
lowest enthalpy phase becomes a cubic Pm-3m structure 
that reaches the remarkable superconducting transition 
temperature of Tc = 139 K at 275 GPa, a Tc much higher 
than any other predicted transition temperature of group 
IV-hydrides. Unfortunately, these results have not been 
experimentally confirmed. 

In this letter, we report on our investigations of the 
disilane system under pressure by using the recently de- 
veloped minima hopping method (MHM) [22, 2^ for the 
prediction of low enthalpy structures. This method has 
been successfully used for global geometry optimization 
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FIG. 2. (Color online) Enthalpy per formula unit of disilane 
as a function of pressure with respect to elements in their 
solid form 2Si(s)+3H2(s). The decomposition enthalpies were 
computed from the predicted structures of hydrogen [sl] and 
high pressure phases of silicon [37, 38]. The disilane structures 
with superscripts ^ and ^ are from this work and from Ref . [2ll, 
respectively. 



in a large variety of applications |2J-|28| . Given only the 
chemical composition of a system, the MHM aims at find- 
ing the global minimum on the enthalpy surface, while 
gradually exploring low lying structures. Moves on the 
enthalpy surface are performed using variable cell shape 
molecular dynamics with initial velocities approximately 
chosen along soft mode directions. The relaxations to lo- 
cal minima are performed by the fast inertia relaxation 
engine (FIRE) [29| taking into account both atomic and 
cell degrees of freedom. 

We performed simulations for cells containing 1,2, and 
3 formula units of disilane Si2H6 under several differ- 
ent pressures between 40-400 GPa. The initial sampling 
of the enthalpy surface was carried out employing the 
MHM together with Lenosky's tight-binding scheme (soj , 
extended to include hydrogen. The most promising 
candidate structures found during the initial sampling 
were further studied 
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|32| at the density functional 
theory (DFT) level using the Perdew-Burke-Erzernhof 
(PBE) exchange-correlation functional ^] and norm- 
conserving HGH-pseudopotentials [s^. The plane- wave 
cut-off energy was set to 1400 eV, and carefully converged 
Monkhorst-Pack k-point meshes [35| with grid spacing 
denser than 2ti x 0.025 A were used. Finally, in order to 
confirm that the tight-binding scheme was able to suffi- 
ciently sample the enthalpy surfaces, we performed MHM 
simulations for selected pressures of 100 GPa, 200 GPa, 
280 GPa, and 320 GPa at the DFT level. 

In Fig. [2] the enthalpy of the different phases found 
in our MHM simulations are shown with respect to 
decomposition towards elemental silicon and hydrogen. 
At pressures above 280 GPa, the Pm-3m phase is fa- 
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FIG. 3. (Color online) Left panel: Phonon band dispersion 
of the Cmcm structure at 200 GPa. Right panel: Calculated 
Eliashberg spectral function a^F(uj) (soUd line) and phonon 
partial density of states (Yellow: silicon, Blue: hydrogen). 



vored, competing with several other structures reported 
by Jin et al [2l| . In addition to the structures reported in 
Ref. [2l| our simulations revealed another low-lying phase 
with P-1 symmetry (Fig. [TJd). However, all these struc- 
tures lie in a very small energy range which is within our 
numerical precision. Furthermore, as seen in Fig.[2l crys- 
talline disilane is enthalpically unstable towards decom- 
position to elemental silicon and hydrogen below 95 GPa. 
A decomposition to silane SiH4 together with elemental 
silicon and hydrogen is enthalpically possible up to pres- 
sures of 190 GPa. This compositional instability could 
pose challenges en route to synthetization of crystalline 
disilane, depending on barrier heights and on the dynam- 
ics of the decomposition. 

Yet another low-enthalpy metallic phase of disilane was 
found during our MHM simulations (see Fig. [T^). It be- 
longs to the Cmcm space group, and is the lowest en- 
thalpy structure up to 280 GPa. At 200 GPa, its con- 
ventional cell parameters are a = 7.965 A, b = 2.705 A, 
and c = 4.728 A, with one silicon atom occupying the 
8e crystallographic site at (0.141,0,0) and three hy- 
drogen atoms occupying 8^, 8^ and 8/ sites at co- 
ordinates (0.293,0.173,0.250), (0.086,0.302,0.250) and 
(0,0.311,0.895), respectively. The hydrogen atoms are 
embedded into a framework of five-fold coordinated sil- 
icon atoms. The average silicon-silicon bond length is 
2.28 A, and each silicon atom is surrounded by six hy- 
drogen atoms at a mean distance of 1.52 A. 

We further characterized the Cmcm structure by 
performing calculations of the phonon-spectrum, the 
electron-phonon coupling and the superconducting tran- 
sition temperature Tc. The phonon spectrum and the 
electron-phonon matrix elements were obtained from 
density- functional perturbation theory [39| . The spectral 
function a'^F{uj) was integrated over the Fermi surface by 
applying the tetrahedron technique. Convergence of the 



TABLE I. Superconducting properties of the Cmcm phase at 
different pressures. The transition temperatures were calcu- 
lated using McMillan's formula. 

Pressure A c<;iog Tc Tc 

^* = 0.1 /i* = 0.13 



100 
140 
160 
200 
220 



0.84 478 

0.68 553 

0.66 556 

0.68 501 

0.76 384 



24.5 
18.0 
16.9 
16.4 
16.4 



20.1 
13.6 
12.8 
12.4 
13.0 



above quantities was ensured byal6xl6xl6 Monkhorst- 
Pack /c-point sampling, and a 4 x 4 x 4 g-point sampling 
for the phonon wave- vectors. The phonon dispersion was 
obtained by Fourier interpolating the computed dynam- 
ical matrices. 

The phonon band dispersion of the Cmcm phase at 
200 GPa can be seen in the left panel of Fig. [3l while the 
the partial phonon density of states is shown in the right 
panel. As expected, the low frequencies (<700cm~^) 
are dominated by the vibrations of the silicon framework 
whereas the high end of the spectrum extending up to 
2300 cm~^ is solely due to the light hydrogen atoms. 
We found the structure to be dynamically stable up to 
220 GPa. However, if the pressure is increased beyond 
225 GPa a dynamical instability arises. The phonon band 
dispersion at 230 GPa can be found in the Supplementary 
Material [40] . It shows an imaginary (plotted as negative) 
frequency at the T-point, indicating an unstable phonon 
mode. Following the eigendisplacements of this mode, 
which are shown in Fig. [It, and then performing a full 
relaxation of the structure leads to another unreported 
stable structure with Cmc2i symmetry. Compared to the 
Cmcm phase, the silicon framework remains essentially 
intact while the hydrogen atoms are slightly displaced, 
partially breaking the symmetry. Due to the strong sim- 
ilarities between the Cmcm and the Cmc2i structures 
we do not expect large differences in their phonons or 
superconducting properties. A similar analysis as above 
has been carried out following a further imaginary fre- 
quency arising at the 6'-point when the pressure is in- 
creased above 260 GPa. The resulting structure found by 
following the corresponding eigendisplacements resulted 
in a structure with Plcl-symmetry (see the Supplemen- 
tary Materials for details on this structure) |4Q| . 

In order to investigate the superconducting proper- 
ties of the Cmcm phase, we use McMillan's approxi- 
mate formula for the superconducting transition tem- 
perature Tc 41, |42[. McMillan's formula requires the 



Eliashberg spectral function a^F(cj), which was ob- 
tained from ab initio calculations performed with the 
abinit code 31, [32[. In the right panel of Fig. [3l the 




FIG. 4. (Color online) Fermi surface of the Cmcm phase at 
several pressures: a) 100 GPa, b) 140 GPa, c) 160 GPa, and 
d) 220 GPa. ^] 



(ii) two intense hydrogen peaks around 1500 cm~^ and 
1600 cm~^, and (iii) high frequency modes of hydrogen 
around 2000 cm~^. 

The superconducting properties of the Cmcm phase 
at several pressures are summarized in Table HI using 
two different values for the Coulomb pseudopotential 
/i* = 0.1 and /i* = 0.13. Assuming the larger of 
those values, the superconducting transition tempera- 
ture Tc is 20.1 K at 100 GPa and decreases to 13.0 K at 
220 GPa. A decreasing Tc with respect to increasing pres- 
sure is a general feature observed in hydrogen rich mate- 



solid lines represent the Eliashberg spectral function of 
the Cmcm phase at 200 GPa. It has three main fea- 
tures: (i) low optical modes of the silicon framework. 



sure IS a general leature observed m nydrogen ricn mate- 
rials [13, US, m • We should emphasize that the Tc of the 
Cmcm phase is approximately one order of magnitude 
smaller than in previously reported structures, and that 
the Cmcm phase is the lowest enthalpy phase. This raises 
serious doubts if high-Tc superconductivity will ever be 
achieved in silane materials under reasonable pressure. 

The superconducting properties of the Cmcm phase 
are strongly linked to its electronic structure. In Fig. SI 
the evolution of the Fermi surface is shown as a func- 
tion of pressure. Three states cross the Fermi surface. 
The first (magenta) and the second (yellow) states cover 
an important portion of the Brillouin zone and remain 
nearly unaltered as the pressure increases, whereas the 
third (cyan) changes substantially. The contribution of 
this third state to the Fermi surface consists of spherical 
regions near the F-point. This band connects two main 
portions of the Fermi surface. Therefore, we can expected 
high superconducting values for low pressures; A = 0.84 
and cjiog = 480 with a Tc of 20 K at 100 GPa. However, 
as the volume of the structure decreases with increas- 
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ing pressure, this sphere-like feature of the Fermi sur- 
face is abruptly reduced. Consequently, at 160 GPa the 
superconducting parameters clearly decrease; A = 0.66, 
'^log = 556 and Tc = 12.4 K. 

In conclusion, we performed a thorough investigation 
of the high pressure phases of disilane using first princi- 
ples calculations. Applying our minima hopping method 
to explore the potential energy surface of disilane, we 
found a metallic structure which is enthalpically favor- 
able compared to the previously proposed structures of 
disilane. Additionally, the systematic study of the su- 
perconducting properties as a function of pressure shows 
that the Cmcm phase possesses a moderate electron- 
phonon coupling, leading to a superconducting transition 
temperature in the 10-20 K range. This result stands in 
sharp contrast with the structures previously proposed of 
disilane under pressure. Moreover, we observed that the 
transition temperature of the Cmcm structure has the 
tendency to decrease monotonically with applied pres- 
sure (for /i* = 0.1), which can be understood by the 
shrinking of a part of the Fermi surface. This decrease 
of the Tc is in agreement with most theoretical and ex- 
perimental results of hydrogen rich materials, including 
silane 17, 43, li^. Certainly, this does not imply that 
superconductivity in hydrogen rich materials is limited 
to relatively low values of Tc, but our results do impose 
strong constraints on the possibility of high- Tc supercon- 
ductors in silicon-hydrogen systems. 
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